/ Shareware Grab Bag / Shareware Grab Bag.iso / 090 / byte0387.arc / DAWSON.ARC / SIMSUBS.C

< prev

next >

MacOS 8.1

Win98

DOS

view JSON data

view as text

+--------+-------------------------+-------------------------+--------+--------+
|00000000| 63 61 74 20 73 69 6d 73 | 75 62 73 2e 63 0d 0a 00 |cat sims|ubs.c...|
|00000010| 00 00 00 2f 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |.../****|********|
|00000020| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000030| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000040| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000050| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 5c 0d 0a 00 00 00 |********|**\.....|
|00000060| 00 2a 20 09 09 09 09 09 | 09 09 09 20 20 20 20 20 |.* .....|...     |
|00000070| 20 20 2a 0d 0a 00 00 00 | 00 2a 20 20 20 20 73 69 |  *.....|.*    si|
|00000080| 6d 73 75 62 73 2e 63 20 | 3d 20 53 75 62 72 6f 75 |msubs.c |= Subrou|
|00000090| 74 69 6e 65 73 20 66 6f | 72 20 73 69 6d 74 65 73 |tines fo|r simtes|
|000000a0| 74 2e 63 20 70 72 6f 67 | 72 61 6d 20 66 6f 72 20 |t.c prog|ram for |
|000000b0| 74 65 73 74 69 6e 67 20 | 53 49 4d 50 50 3a 20 20 |testing |SIMPP:  |
|000000c0| 2a 0d 0a 00 00 00 00 2a | 20 20 20 20 46 6f 72 20 |*......*|    For |
|000000d0| 53 69 6d 70 6c 65 20 49 | 4d 61 67 65 20 50 72 6f |Simple I|Mage Pro|
|000000e0| 63 65 73 73 69 6e 67 20 | 50 61 63 6b 61 67 65 2e |cessing |Package.|
|000000f0| 09 09 09 20 20 20 20 20 | 20 20 2a 0d 0a 00 00 00 |...     |  *.....|
|00000100| 00 2a 20 20 20 20 43 6f | 70 79 72 69 67 68 74 20 |.*    Co|pyright |
|00000110| 28 63 29 20 31 39 38 37 | 2c 20 42 65 6e 6a 61 6d |(c) 1987|, Benjam|
|00000120| 69 6e 20 4d 2e 20 44 61 | 77 73 6f 6e 2e 09 09 09 |in M. Da|wson....|
|00000130| 20 20 20 20 20 20 20 2a | 0d 0a 00 00 00 00 2a 20 |       *|......* |
|00000140| 20 20 20 20 20 45 64 69 | 74 20 56 65 72 73 69 6f |     Edi|t Versio|
|00000150| 6e 3a 20 31 2e 31 20 3a | 20 4a 61 6e 2d 32 39 2d |n: 1.1 :| Jan-29-|
|00000160| 38 37 09 09 09 09 20 20 | 20 20 20 20 20 2a 0d 0a |87....  |     *..|
|00000170| 00 00 00 00 2a 20 09 09 | 09 09 09 09 09 09 20 20 |....* ..|......  |
|00000180| 20 20 20 20 20 2a 0d 0a | 00 00 00 00 5c 2a 2a 2a |     *..|....\***|
|00000190| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000001a0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000001b0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000001c0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000001d0| 2a 2a 2a 2f 0d 0a 00 00 | 00 00 0d 0a 00 00 00 00 |***/....|........|
|000001e0| 23 69 6e 63 6c 75 64 65 | 20 3c 73 74 64 69 6f 2e |#include| <stdio.|
|000001f0| 68 3e 0d 0a 00 00 00 00 | 23 69 6e 63 6c 75 64 65 |h>......|#include|
|00000200| 20 22 73 69 6d 70 70 2e | 68 22 0d 0a 00 00 00 00 | "simpp.|h"......|
|00000210| 0d 0a 00 00 00 00 65 78 | 74 65 72 6e 20 64 6f 75 |......ex|tern dou|
|00000220| 62 6c 65 20 65 78 70 28 | 29 3b 09 09 09 2f 2a 20 |ble exp(|);.../* |
|00000230| 52 45 51 55 49 52 45 53 | 20 4d 41 54 48 20 4c 49 |REQUIRES| MATH LI|
|00000240| 42 52 41 52 59 21 20 2a | 2f 0d 0a 00 00 00 00 73 |BRARY! *|/......s|
|00000250| 74 61 74 69 63 20 50 49 | 58 45 4c 20 74 72 61 6e |tatic PI|XEL tran|
|00000260| 5b 50 49 58 45 4c 5f 53 | 49 5a 45 5d 20 3d 20 30 |[PIXEL_S|IZE] = 0|
|00000270| 3b 0d 0a 00 00 00 00 0d | 0a 00 00 00 00 2f 2a 20 |;.......|...../* |
|00000280| 3d 3d 3d 3d 3d 3d 3d 3d | 3d 3d 3d 3d 3d 3d 3d 3d |========|========|
|00000290| 3d 3d 3d 3d 20 41 64 64 | 69 74 69 6f 6e 61 6c 20 |==== Add|itional |
|000002a0| 50 6f 69 6e 74 20 72 6f | 75 74 69 6e 65 73 20 3d |Point ro|utines =|
|000002b0| 3d 3d 3d 3d 3d 3d 3d 3d | 3d 3d 3d 3d 3d 3d 3d 3d |========|========|
|000002c0| 3d 3d 20 2a 2f 0d 0a 00 | 00 00 00 0d 0a 00 00 00 |== */...|........|
|000002d0| 00 2f 2a 20 2b 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b |./* ++++|++++++++|
|000002e0| 2b 2b 2b 2b 2b 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b |++++++++|++++++++|
|000002f0| 2b 2b 2b 2b 2b 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b |++++++++|++++++++|
|00000300| 2b 2b 2b 2b 2b 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b |++++++++|++++++++|
|00000310| 2b 2b 2b 2b 20 2a 2f 0d | 0a 00 00 00 00 56 4f 49 |++++ */.|.....VOI|
|00000320| 44 20 6e 65 67 61 74 65 | 28 78 2c 79 2c 64 78 2c |D negate|(x,y,dx,|
|00000330| 64 79 29 09 09 2f 2a 20 | 4e 65 67 61 74 65 20 61 |dy)../* |Negate a|
|00000340| 6e 20 61 72 65 61 20 2a | 2f 0d 0a 00 00 00 00 69 |n area *|/......i|
|00000350| 6e 74 20 78 2c 79 2c 64 | 78 2c 64 79 3b 09 09 09 |nt x,y,d|x,dy;...|
|00000360| 2f 2a 20 41 72 65 61 20 | 74 6f 20 6e 65 67 61 74 |/* Area |to negat|
|00000370| 65 20 2a 2f 0d 0a 00 00 | 00 00 7b 0d 0a 00 00 00 |e */....|..{.....|
|00000380| 00 09 72 65 67 69 73 74 | 65 72 20 69 6e 74 20 69 |..regist|er int i|
|00000390| 3b 0d 0a 00 00 00 00 0d | 0a 00 00 00 00 09 66 6f |;.......|......fo|
|000003a0| 72 20 28 69 20 3d 20 30 | 20 3b 20 69 20 3c 20 50 |r (i = 0| ; i < P|
|000003b0| 49 58 45 4c 5f 53 49 5a | 45 20 3b 20 69 2b 2b 29 |IXEL_SIZ|E ; i++)|
|000003c0| 09 2f 2a 20 53 65 74 20 | 75 70 20 74 72 61 6e 73 |./* Set |up trans|
|000003d0| 66 6f 72 6d 61 74 69 6f | 6e 20 2a 2f 0d 0a 00 00 |formatio|n */....|
|000003e0| 00 00 09 09 74 72 61 6e | 5b 69 5d 20 3d 20 4d 41 |....tran|[i] = MA|
|000003f0| 58 50 49 58 20 2d 20 28 | 50 49 58 45 4c 29 69 3b |XPIX - (|PIXEL)i;|
|00000400| 0d 0a 00 00 00 00 0d 0a | 00 00 00 00 09 70 74 72 |........|.....ptr|
|00000410| 61 6e 73 66 6f 72 6d 28 | 78 2c 79 2c 64 78 2c 64 |ansform(|x,y,dx,d|
|00000420| 79 2c 74 72 61 6e 29 3b | 09 09 2f 2a 20 41 70 70 |y,tran);|../* App|
|00000430| 6c 79 20 69 74 20 2a 2f | 0d 0a 00 00 00 00 7d 0d |ly it */|......}.|
|00000440| 0a 00 00 00 00 0d 0a 00 | 00 00 00 2f 2a 20 2b 2b |........|.../* ++|
|00000450| 2b 2b 2b 2b 2b 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b |++++++++|++++++++|
|00000460| 2b 2b 2b 2b 2b 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b |++++++++|++++++++|
|00000470| 2b 2b 2b 2b 2b 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b |++++++++|++++++++|
|00000480| 2b 2b 2b 2b 2b 2b 2b 2b | 2b 2b 2b 2b 2b 2b 20 2a |++++++++|++++++ *|
|00000490| 2f 09 0d 0a 00 00 00 00 | 56 4f 49 44 20 62 72 69 |/.......|VOID bri|
|000004a0| 67 68 74 65 6e 28 78 2c | 79 2c 64 78 2c 64 79 2c |ghten(x,|y,dx,dy,|
|000004b0| 76 61 6c 29 09 2f 2a 20 | 62 72 69 67 68 74 65 6e |val)./* |brighten|
|000004c0| 20 61 6e 20 61 72 65 61 | 2e 20 20 44 69 6d 20 69 | an area|.  Dim i|
|000004d0| 66 20 76 61 6c 20 3c 20 | 30 20 2a 2f 0d 0a 00 00 |f val < |0 */....|
|000004e0| 00 00 69 6e 74 20 78 2c | 79 2c 64 78 2c 64 79 3b |..int x,|y,dx,dy;|
|000004f0| 09 09 09 2f 2a 20 61 72 | 65 61 20 74 6f 20 62 72 |.../* ar|ea to br|
|00000500| 69 67 68 74 65 6e 20 2a | 2f 0d 0a 00 00 00 00 69 |ighten *|/......i|
|00000510| 6e 74 20 76 61 6c 3b 09 | 09 09 2f 2a 20 56 61 6c |nt val;.|../* Val|
|00000520| 75 65 20 74 6f 20 61 64 | 64 20 74 6f 20 61 72 65 |ue to ad|d to are|
|00000530| 61 20 2a 2f 0d 0a 00 00 | 00 00 7b 0d 0a 00 00 00 |a */....|..{.....|
|00000540| 00 09 72 65 67 69 73 74 | 65 72 20 69 6e 74 20 69 |..regist|er int i|
|00000550| 3b 0d 0a 00 00 00 00 0d | 0a 00 00 00 00 09 69 66 |;.......|......if|
|00000560| 20 28 76 61 6c 20 3e 20 | 30 29 20 7b 09 2f 2a 20 | (val > |0) {./* |
|00000570| 42 52 49 47 48 54 45 4e | 21 21 20 2a 2f 0d 0a 00 |BRIGHTEN|!! */...|
|00000580| 00 00 00 09 20 20 66 6f | 72 20 28 69 20 3d 20 30 |....  fo|r (i = 0|
|00000590| 20 3b 20 69 20 3c 20 50 | 49 58 45 4c 5f 53 49 5a | ; i < P|IXEL_SIZ|
|000005a0| 45 20 3b 20 69 2b 2b 29 | 20 7b 09 2f 2a 20 53 65 |E ; i++)| {./* Se|
|000005b0| 74 20 75 70 20 74 72 61 | 6e 73 66 6f 72 6d 20 2a |t up tra|nsform *|
|000005c0| 2f 0d 0a 00 00 00 00 09 | 09 69 66 20 28 28 69 2b |/.......|.if ((i+|
|000005d0| 76 61 6c 29 20 3e 20 28 | 75 6e 73 69 67 6e 65 64 |val) > (|unsigned|
|000005e0| 20 69 6e 74 29 4d 41 58 | 50 49 58 29 20 74 72 61 | int)MAX|PIX) tra|
|000005f0| 6e 5b 69 5d 20 3d 20 4d | 41 58 50 49 58 3b 0d 0a |n[i] = M|AXPIX;..|
|00000600| 00 00 00 00 09 09 65 6c | 73 65 20 74 72 61 6e 5b |......el|se tran[|
|00000610| 69 5d 20 3d 20 28 50 49 | 58 45 4c 29 28 69 2b 76 |i] = (PI|XEL)(i+v|
|00000620| 61 6c 29 3b 0d 0a 00 00 | 00 00 09 20 20 7d 0d 0a |al);....|...  }..|
|00000630| 00 00 00 00 09 7d 0d 0a | 00 00 00 00 09 65 6c 73 |.....}..|.....els|
|00000640| 65 20 7b 09 09 2f 2a 20 | 44 49 4d 20 21 21 20 2a |e {../* |DIM !! *|
|00000650| 2f 0d 0a 00 00 00 00 09 | 20 20 66 6f 72 20 28 69 |/.......|  for (i|
|00000660| 20 3d 20 30 20 3b 20 69 | 20 3c 20 50 49 58 45 4c | = 0 ; i| < PIXEL|
|00000670| 5f 53 49 5a 45 20 3b 20 | 69 2b 2b 29 20 7b 09 2f |_SIZE ; |i++) {./|
|00000680| 2a 20 53 65 74 20 75 70 | 20 74 72 61 6e 73 66 6f |* Set up| transfo|
|00000690| 72 6d 20 2a 2f 0d 0a 00 | 00 00 00 09 09 69 66 20 |rm */...|.....if |
|000006a0| 28 28 69 2b 76 61 6c 29 | 20 3c 20 28 69 6e 74 29 |((i+val)| < (int)|
|000006b0| 4d 49 4e 50 49 58 29 20 | 74 72 61 6e 5b 69 5d 20 |MINPIX) |tran[i] |
|000006c0| 3d 20 4d 49 4e 50 49 58 | 3b 0d 0a 00 00 00 00 09 |= MINPIX|;.......|
|000006d0| 09 65 6c 73 65 20 74 72 | 61 6e 5b 69 5d 20 3d 20 |.else tr|an[i] = |
|000006e0| 28 50 49 58 45 4c 29 28 | 69 2b 76 61 6c 29 3b 0d |(PIXEL)(|i+val);.|
|000006f0| 0a 00 00 00 00 09 20 20 | 7d 0d 0a 00 00 00 00 09 |......  |}.......|
|00000700| 7d 0d 0a 00 00 00 00 0d | 0a 00 00 00 00 09 70 74 |}.......|......pt|
|00000710| 72 61 6e 73 66 6f 72 6d | 28 78 2c 79 2c 64 78 2c |ransform|(x,y,dx,|
|00000720| 64 79 2c 74 72 61 6e 29 | 3b 09 09 2f 2a 20 41 70 |dy,tran)|;../* Ap|
|00000730| 70 6c 79 20 69 74 20 2a | 2f 0d 0a 00 00 00 00 7d |ply it *|/......}|
|00000740| 0d 0a 00 00 00 00 0d 0a | 00 00 00 00 2f 2a 20 2b |........|..../* +|
|00000750| 2b 2b 2b 2b 2b 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b |++++++++|++++++++|
|00000760| 2b 2b 2b 2b 2b 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b |++++++++|++++++++|
|00000770| 2b 2b 2b 2b 2b 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b |++++++++|++++++++|
|00000780| 2b 2b 2b 2b 2b 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 20 |++++++++|+++++++ |
|00000790| 2a 2f 0d 0a 00 00 00 00 | 2f 2a 20 22 42 75 72 6e |*/......|/* "Burn|
|000007a0| 69 6e 67 22 20 69 73 20 | 61 20 70 68 6f 74 6f 67 |ing" is |a photog|
|000007b0| 72 61 70 68 69 63 20 74 | 65 63 68 6e 69 71 75 65 |raphic t|echnique|
|000007c0| 20 74 6f 20 65 6e 68 61 | 6e 63 65 20 74 68 65 20 | to enha|nce the |
|000007d0| 65 78 70 6f 73 75 72 65 | 20 61 6e 64 0d 0a 00 00 |exposure| and....|
|000007e0| 00 00 20 2a 20 68 65 6e | 63 65 20 74 68 65 20 63 |.. * hen|ce the c|
|000007f0| 6f 6e 74 72 61 73 74 20 | 69 6e 20 73 65 6c 65 63 |ontrast |in selec|
|00000800| 74 69 76 65 20 69 6d 61 | 67 65 20 61 72 65 61 73 |tive ima|ge areas|
|00000810| 2e 20 20 54 68 69 73 20 | 69 73 20 64 6f 6e 65 20 |.  This |is done |
|00000820| 62 79 20 6d 6f 76 69 6e | 67 0d 0a 00 00 00 00 20 |by movin|g...... |
|00000830| 2a 20 61 20 6d 61 73 6b | 20 77 69 74 68 20 61 20 |* a mask| with a |
|00000840| 63 75 74 6f 75 74 20 28 | 75 73 75 61 6c 6c 79 20 |cutout (|usually |
|00000850| 61 20 63 69 72 63 6c 65 | 29 20 61 72 6f 75 6e 64 |a circle|) around|
|00000860| 20 74 68 65 20 61 72 65 | 61 20 77 68 69 6c 65 20 | the are|a while |
|00000870| 69 74 0d 0a 00 00 00 00 | 20 2a 20 69 73 20 62 65 |it......| * is be|
|00000880| 69 6e 67 20 65 6e 6c 61 | 72 67 65 64 2e 20 20 57 |ing enla|rged.  W|
|00000890| 65 20 61 70 70 72 6f 78 | 69 6d 61 74 65 20 74 68 |e approx|imate th|
|000008a0| 69 73 20 70 72 6f 63 65 | 64 75 72 65 20 75 73 69 |is proce|dure usi|
|000008b0| 6e 67 20 61 20 67 75 61 | 73 73 69 61 6e 0d 0a 00 |ng a gua|ssian...|
|000008c0| 00 00 00 20 2a 20 73 68 | 61 70 65 64 20 65 6e 68 |... * sh|aped enh|
|000008d0| 61 6e 63 65 6d 65 6e 74 | 20 6f 66 20 74 68 65 20 |ancement| of the |
|000008e0| 63 6f 6e 74 72 61 73 74 | 2e 0d 0a 00 00 00 00 20 |contrast|....... |
|000008f0| 2a 2f 0d 0a 00 00 00 00 | 56 4f 49 44 20 67 61 75 |*/......|VOID gau|
|00000900| 73 73 5f 62 75 72 6e 28 | 78 2c 79 2c 64 78 2c 64 |ss_burn(|x,y,dx,d|
|00000910| 79 2c 6b 2c 6c 2c 6d 29 | 09 2f 2a 20 47 61 75 73 |y,k,l,m)|./* Gaus|
|00000920| 73 69 61 6e 20 22 42 75 | 72 6e 22 20 61 6e 20 61 |sian "Bu|rn" an a|
|00000930| 72 65 61 20 2a 2f 0d 0a | 00 00 00 00 69 6e 74 20 |rea */..|....int |
|00000940| 78 2c 79 2c 64 78 2c 64 | 79 3b 09 09 2f 2a 20 41 |x,y,dx,d|y;../* A|
|00000950| 72 65 61 20 74 6f 20 62 | 75 72 6e 20 2a 2f 0d 0a |rea to b|urn */..|
|00000960| 00 00 00 00 64 6f 75 62 | 6c 65 20 6b 3b 09 09 2f |....doub|le k;../|
|00000970| 2a 20 45 78 70 61 6e 73 | 69 6f 6e 20 28 62 75 72 |* Expans|ion (bur|
|00000980| 6e 29 20 66 61 63 74 6f | 72 20 2a 2f 0d 0a 00 00 |n) facto|r */....|
|00000990| 00 00 64 6f 75 62 6c 65 | 20 6c 3b 09 09 2f 2a 20 |..double| l;../* |
|000009a0| 58 20 61 6e 64 20 59 20 | 73 70 61 63 65 20 63 6f |X and Y |space co|
|000009b0| 6e 73 74 61 6e 74 73 20 | 2a 2f 0d 0a 00 00 00 00 |nstants |*/......|
|000009c0| 64 6f 75 62 6c 65 20 6d | 3b 09 09 2f 2a 20 4f 66 |double m|;../* Of|
|000009d0| 66 73 65 74 20 66 61 63 | 74 6f 72 20 2a 2f 0d 0a |fset fac|tor */..|
|000009e0| 00 00 00 00 7b 0d 0a 00 | 00 00 00 09 72 65 67 69 |....{...|....regi|
|000009f0| 73 74 65 72 20 69 6e 74 | 20 69 3b 0d 0a 00 00 00 |ster int| i;.....|
|00000a00| 00 09 64 6f 75 62 6c 65 | 20 78 78 2c 20 79 79 2c |..double| xx, yy,|
|00000a10| 20 7a 3b 0d 0a 00 00 00 | 00 0d 0a 00 00 00 00 09 | z;.....|........|
|00000a20| 79 79 20 3d 20 2d 28 64 | 6f 75 62 6c 65 29 64 79 |yy = -(d|ouble)dy|
|00000a30| 2f 32 2e 30 3b 0d 0a 00 | 00 00 00 0d 0a 00 00 00 |/2.0;...|........|
|00000a40| 00 09 77 68 69 6c 65 20 | 28 64 79 2d 2d 29 20 7b |..while |(dy--) {|
|00000a50| 0d 0a 00 00 00 00 09 20 | 20 78 78 20 3d 20 2d 28 |....... | xx = -(|
|00000a60| 64 6f 75 62 6c 65 29 64 | 78 2f 32 2e 30 3b 0d 0a |double)d|x/2.0;..|
|00000a70| 00 00 00 00 09 20 20 66 | 6f 72 20 28 69 20 3d 20 |.....  f|or (i = |
|00000a80| 78 20 3b 20 69 20 3c 20 | 78 2b 64 78 20 3b 20 69 |x ; i < |x+dx ; i|
|00000a90| 2b 2b 29 20 7b 0d 0a 00 | 00 00 00 09 09 7a 20 3d |++) {...|.....z =|
|00000aa0| 20 28 64 6f 75 62 6c 65 | 29 72 65 61 64 5f 70 69 | (double|)read_pi|
|00000ab0| 78 65 6c 28 69 2c 79 29 | 2a 6b 2a 65 78 70 28 2d |xel(i,y)|*k*exp(-|
|00000ac0| 28 28 78 78 2a 78 78 20 | 2b 20 79 79 2a 79 79 29 |((xx*xx |+ yy*yy)|
|00000ad0| 2f 6c 29 29 20 2d 20 6d | 3b 0d 0a 00 00 00 00 09 |/l)) - m|;.......|
|00000ae0| 09 69 66 20 28 7a 20 3c | 20 30 2e 30 29 20 77 72 |.if (z <| 0.0) wr|
|00000af0| 69 74 65 5f 70 69 78 65 | 6c 28 69 2c 79 2c 4d 49 |ite_pixe|l(i,y,MI|
|00000b00| 4e 50 49 58 29 3b 0d 0a | 00 00 00 00 09 09 65 6c |NPIX);..|......el|
|00000b10| 73 65 20 69 66 20 28 7a | 20 3e 20 28 64 6f 75 62 |se if (z| > (doub|
|00000b20| 6c 65 29 4d 41 58 50 49 | 58 29 20 77 72 69 74 65 |le)MAXPI|X) write|
|00000b30| 5f 70 69 78 65 6c 28 69 | 2c 79 2c 4d 41 58 50 49 |_pixel(i|,y,MAXPI|
|00000b40| 58 29 3b 0d 0a 00 00 00 | 00 09 09 65 6c 73 65 20 |X);.....|...else |
|00000b50| 77 72 69 74 65 5f 70 69 | 78 65 6c 28 69 2c 79 2c |write_pi|xel(i,y,|
|00000b60| 20 28 50 49 58 45 4c 29 | 7a 29 3b 0d 0a 00 00 00 | (PIXEL)|z);.....|
|00000b70| 00 09 09 78 78 20 2b 3d | 20 31 2e 30 3b 0d 0a 00 |...xx +=| 1.0;...|
|00000b80| 00 00 00 09 20 20 7d 0d | 0a 00 00 00 00 09 20 20 |....  }.|......  |
|00000b90| 79 79 20 2b 3d 20 31 2e | 30 3b 0d 0a 00 00 00 00 |yy += 1.|0;......|
|00000ba0| 09 20 20 79 2b 2b 3b 0d | 0a 00 00 00 00 09 7d 0d |.  y++;.|......}.|
|00000bb0| 0a 00 00 00 00 7d 0d 0a | 00 00 00 00 0d 0a 00 00 |.....}..|........|
|00000bc0| 00 00 2f 2a 20 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b |../* +++|++++++++|
|00000bd0| 2b 2b 2b 2b 2b 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b |++++++++|++++++++|
|00000be0| 2b 2b 2b 2b 2b 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b |++++++++|++++++++|
|00000bf0| 2b 2b 2b 2b 2b 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b |++++++++|++++++++|
|00000c00| 2b 2b 2b 2b 2b 20 2a 2f | 0d 0a 00 00 00 00 56 4f |+++++ */|......VO|
|00000c10| 49 44 20 62 69 6e 61 72 | 69 7a 65 28 78 2c 79 2c |ID binar|ize(x,y,|
|00000c20| 64 78 2c 64 79 2c 74 68 | 72 65 73 68 6f 6c 64 29 |dx,dy,th|reshold)|
|00000c30| 09 09 2f 2a 20 42 69 6e | 61 72 69 7a 65 20 61 72 |../* Bin|arize ar|
|00000c40| 65 61 20 2a 2f 0d 0a 00 | 00 00 00 69 6e 74 20 78 |ea */...|...int x|
|00000c50| 2c 79 2c 64 78 2c 64 79 | 3b 09 09 2f 2a 20 41 72 |,y,dx,dy|;../* Ar|
|00000c60| 65 61 20 74 6f 20 62 69 | 6e 61 72 69 7a 65 20 2a |ea to bi|narize *|
|00000c70| 2f 0d 0a 00 00 00 00 69 | 6e 74 20 74 68 72 65 73 |/......i|nt thres|
|00000c80| 68 6f 6c 64 3b 09 09 2f | 2a 20 54 68 72 65 73 68 |hold;../|* Thresh|
|00000c90| 6f 6c 64 20 76 61 6c 75 | 65 20 2a 2f 0d 0a 00 00 |old valu|e */....|
|00000ca0| 00 00 7b 0d 0a 00 00 00 | 00 09 72 65 67 69 73 74 |..{.....|..regist|
|00000cb0| 65 72 20 69 6e 74 20 69 | 3b 0d 0a 00 00 00 00 0d |er int i|;.......|
|00000cc0| 0a 00 00 00 00 09 66 6f | 72 20 28 69 20 3d 20 30 |......fo|r (i = 0|
|00000cd0| 20 3b 20 69 20 3c 3d 20 | 74 68 72 65 73 68 6f 6c | ; i <= |threshol|
|00000ce0| 64 20 3b 20 69 2b 2b 29 | 20 74 72 61 6e 5b 69 5d |d ; i++)| tran[i]|
|00000cf0| 20 3d 20 4d 49 4e 50 49 | 58 3b 0d 0a 00 00 00 00 | = MINPI|X;......|
|00000d00| 09 66 6f 72 20 28 69 20 | 3d 20 74 68 72 65 73 68 |.for (i |= thresh|
|00000d10| 6f 6c 64 20 3b 20 69 20 | 3c 20 50 49 58 45 4c 5f |old ; i |< PIXEL_|
|00000d20| 53 49 5a 45 20 3b 20 69 | 2b 2b 29 20 74 72 61 6e |SIZE ; i|++) tran|
|00000d30| 5b 69 5d 20 3d 20 4d 41 | 58 50 49 58 3b 0d 0a 00 |[i] = MA|XPIX;...|
|00000d40| 00 00 00 09 70 74 72 61 | 6e 73 66 6f 72 6d 28 78 |....ptra|nsform(x|
|00000d50| 2c 79 2c 64 78 2c 64 79 | 2c 74 72 61 6e 29 3b 0d |,y,dx,dy|,tran);.|
|00000d60| 0a 00 00 00 00 7d 0d 0a | 00 00 00 00 0d 0a 00 00 |.....}..|........|
|00000d70| 00 00 2f 2a 20 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b |../* +++|++++++++|
|00000d80| 2b 2b 2b 2b 2b 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b |++++++++|++++++++|
|00000d90| 2b 2b 2b 2b 2b 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b |++++++++|++++++++|
|00000da0| 2b 2b 2b 2b 2b 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b |++++++++|++++++++|
|00000db0| 2b 2b 2b 2b 2b 20 2a 2f | 0d 0a 00 00 00 00 56 4f |+++++ */|......VO|
|00000dc0| 49 44 20 63 73 74 72 65 | 74 63 68 28 78 2c 79 2c |ID cstre|tch(x,y,|
|00000dd0| 64 78 2c 64 79 2c 63 6c | 69 70 29 09 09 2f 2a 20 |dx,dy,cl|ip)../* |
|00000de0| 53 74 72 65 74 63 68 20 | 63 6f 6e 74 72 61 73 74 |Stretch |contrast|
|00000df0| 20 69 6e 20 61 72 65 61 | 20 2a 2f 0d 0a 00 00 00 | in area| */.....|
|00000e00| 00 69 6e 74 20 78 2c 79 | 2c 64 78 2c 64 79 3b 09 |.int x,y|,dx,dy;.|
|00000e10| 09 09 09 2f 2a 20 41 72 | 65 61 20 74 6f 20 65 6e |.../* Ar|ea to en|
|00000e20| 68 61 6e 63 65 20 2a 2f | 0d 0a 00 00 00 00 69 6e |hance */|......in|
|00000e30| 74 20 63 6c 69 70 3b 09 | 09 09 09 2f 2a 20 4d 69 |t clip;.|.../* Mi|
|00000e40| 6e 69 6d 75 6d 20 62 69 | 6e 20 63 6f 75 6e 74 20 |nimum bi|n count |
|00000e50| 2a 2f 0d 0a 00 00 00 00 | 7b 0d 0a 00 00 00 00 09 |*/......|{.......|
|00000e60| 72 65 67 69 73 74 65 72 | 20 69 6e 74 20 69 3b 0d |register| int i;.|
|00000e70| 0a 00 00 00 00 09 69 6e | 74 20 6c 6f 77 5f 62 69 |......in|t low_bi|
|00000e80| 6e 2c 68 69 67 68 5f 62 | 69 6e 3b 09 09 2f 2a 20 |n,high_b|in;../* |
|00000e90| 48 69 73 74 6f 67 72 61 | 6d 20 68 69 67 68 2c 20 |Histogra|m high, |
|00000ea0| 6c 6f 77 20 62 69 6e 73 | 20 2a 2f 0d 0a 00 00 00 |low bins| */.....|
|00000eb0| 00 09 64 6f 75 62 6c 65 | 20 73 74 65 70 2c 73 74 |..double| step,st|
|00000ec0| 65 70 5f 76 61 6c 75 65 | 3b 09 09 2f 2a 20 46 6f |ep_value|;../* Fo|
|00000ed0| 72 20 68 69 73 74 6f 67 | 72 61 6d 20 65 71 75 61 |r histog|ram equa|
|00000ee0| 6c 69 7a 61 74 69 6f 6e | 20 2a 2f 0d 0a 00 00 00 |lization| */.....|
|00000ef0| 00 09 6c 6f 6e 67 20 68 | 5b 50 49 58 45 4c 5f 53 |..long h|[PIXEL_S|
|00000f00| 49 5a 45 5d 3b 09 09 2f | 2a 20 48 69 73 74 6f 67 |IZE];../|* Histog|
|00000f10| 72 61 6d 20 61 72 72 61 | 79 20 2a 2f 0d 0a 00 00 |ram arra|y */....|
|00000f20| 00 00 0d 0a 00 00 00 00 | 09 68 69 73 74 6f 67 72 |........|.histogr|
|00000f30| 61 6d 28 78 2c 79 2c 64 | 78 2c 64 79 2c 68 29 3b |am(x,y,d|x,dy,h);|
|00000f40| 09 09 2f 2a 20 48 69 73 | 74 6f 67 72 61 6d 20 61 |../* His|togram a|
|00000f50| 72 65 61 20 2a 2f 0d 0a | 00 00 00 00 0d 0a 00 00 |rea */..|........|
|00000f60| 00 00 09 63 6c 69 70 5f | 68 69 73 74 6f 28 68 2c |...clip_|histo(h,|
|00000f70| 63 6c 69 70 2c 26 6c 6f | 77 5f 62 69 6e 2c 26 68 |clip,&lo|w_bin,&h|
|00000f80| 69 67 68 5f 62 69 6e 29 | 3b 09 2f 2a 20 67 65 74 |igh_bin)|;./* get|
|00000f90| 20 6d 69 6e 2c 6d 61 78 | 20 68 69 73 74 6f 20 62 | min,max| histo b|
|00000fa0| 69 6e 73 20 2a 2f 0d 0a | 00 00 00 00 09 70 72 69 |ins */..|.....pri|
|00000fb0| 6e 74 66 28 22 63 73 74 | 72 65 74 63 68 3a 20 4d |ntf("cst|retch: M|
|00000fc0| 69 6e 75 6d 75 6d 20 62 | 69 6e 20 3d 20 25 64 2c |inumum b|in = %d,|
|00000fd0| 20 6d 61 78 69 6d 75 6d | 20 3d 20 25 64 5c 6e 22 | maximum| = %d\n"|
|00000fe0| 2c 6c 6f 77 5f 62 69 6e | 2c 68 69 67 68 5f 62 69 |,low_bin|,high_bi|
|00000ff0| 6e 29 3b 0d 0a 00 00 00 | 00 0d 0a 00 00 00 00 09 |n);.....|........|
|00001000| 73 74 65 70 20 3d 20 28 | 64 6f 75 62 6c 65 29 50 |step = (|double)P|
|00001010| 49 58 45 4c 5f 53 49 5a | 45 2f 28 64 6f 75 62 6c |IXEL_SIZ|E/(doubl|
|00001020| 65 29 28 68 69 67 68 5f | 62 69 6e 2d 6c 6f 77 5f |e)(high_|bin-low_|
|00001030| 62 69 6e 2b 31 29 3b 2f | 2a 20 73 74 65 70 20 73 |bin+1);/|* step s|
|00001040| 69 7a 65 20 2a 2f 0d 0a | 00 00 00 00 09 73 74 65 |ize */..|.....ste|
|00001050| 70 5f 76 61 6c 75 65 20 | 3d 20 30 2e 30 3b 09 09 |p_value |= 0.0;..|
|00001060| 2f 2a 20 53 74 65 70 20 | 76 61 6c 75 65 20 2a 2f |/* Step |value */|
|00001070| 0d 0a 00 00 00 00 2f 2a | 20 56 61 6c 75 65 73 20 |....../*| Values |
|00001080| 62 65 6c 6f 77 20 6c 6f | 77 5f 62 69 6e 20 61 72 |below lo|w_bin ar|
|00001090| 65 20 73 65 74 20 74 6f | 20 6d 69 6e 69 6d 75 6d |e set to| minimum|
|000010a0| 20 70 69 78 65 6c 20 76 | 61 6c 75 65 20 2a 2f 0d | pixel v|alue */.|
|000010b0| 0a 00 00 00 00 09 66 6f | 72 20 28 69 20 3d 20 30 |......fo|r (i = 0|
|000010c0| 20 3b 20 69 20 3c 20 6c | 6f 77 5f 62 69 6e 20 3b | ; i < l|ow_bin ;|
|000010d0| 20 69 2b 2b 29 20 74 72 | 61 6e 5b 69 5d 20 3d 20 | i++) tr|an[i] = |
|000010e0| 4d 49 4e 50 49 58 3b 0d | 0a 00 00 00 00 2f 2a 20 |MINPIX;.|...../* |
|000010f0| 56 61 6c 75 65 73 20 62 | 65 74 77 65 65 6e 20 6c |Values b|etween l|
|00001100| 6f 77 5f 62 69 6e 20 61 | 6e 64 20 68 69 67 68 5f |ow_bin a|nd high_|
|00001110| 62 69 6e 20 61 72 65 20 | 73 74 72 65 74 63 68 65 |bin are |stretche|
|00001120| 64 20 74 6f 20 72 61 6e | 67 65 20 66 72 6f 6d 20 |d to ran|ge from |
|00001130| 4d 49 4e 50 49 58 20 74 | 6f 0d 0a 00 00 00 00 20 |MINPIX t|o...... |
|00001140| 20 20 4d 41 58 50 49 58 | 20 2a 2f 0d 0a 00 00 00 |  MAXPIX| */.....|
|00001150| 00 09 66 6f 72 20 28 69 | 20 3d 20 6c 6f 77 5f 62 |..for (i| = low_b|
|00001160| 69 6e 20 3b 20 69 20 3c | 3d 20 68 69 67 68 5f 62 |in ; i <|= high_b|
|00001170| 69 6e 20 3b 20 69 2b 2b | 29 20 7b 0d 0a 00 00 00 |in ; i++|) {.....|
|00001180| 00 09 09 74 72 61 6e 5b | 69 5d 20 3d 20 28 50 49 |...tran[|i] = (PI|
|00001190| 58 45 4c 29 73 74 65 70 | 5f 76 61 6c 75 65 3b 0d |XEL)step|_value;.|
|000011a0| 0a 00 00 00 00 09 09 73 | 74 65 70 5f 76 61 6c 75 |.......s|tep_valu|
|000011b0| 65 20 2b 3d 20 73 74 65 | 70 3b 0d 0a 00 00 00 00 |e += ste|p;......|
|000011c0| 09 7d 0d 0a 00 00 00 00 | 2f 2a 20 56 61 6c 75 65 |.}......|/* Value|
|000011d0| 73 20 61 62 6f 76 65 20 | 68 69 67 68 5f 62 69 6e |s above |high_bin|
|000011e0| 20 61 72 65 20 73 65 74 | 20 74 6f 20 4d 41 58 50 | are set| to MAXP|
|000011f0| 49 58 20 2a 2f 0d 0a 00 | 00 00 00 09 66 6f 72 20 |IX */...|....for |
|00001200| 28 69 20 3d 20 68 69 67 | 68 5f 62 69 6e 2b 31 20 |(i = hig|h_bin+1 |
|00001210| 3b 20 69 20 3c 20 50 49 | 58 45 4c 5f 53 49 5a 45 |; i < PI|XEL_SIZE|
|00001220| 20 3b 20 69 2b 2b 29 20 | 74 72 61 6e 5b 69 5d 20 | ; i++) |tran[i] |
|00001230| 3d 20 4d 41 58 50 49 58 | 3b 0d 0a 00 00 00 00 0d |= MAXPIX|;.......|
|00001240| 0a 00 00 00 00 2f 2a 20 | 54 72 61 6e 73 66 6f 72 |...../* |Transfor|
|00001250| 6d 20 61 72 65 61 20 62 | 79 20 73 74 72 65 74 63 |m area b|y stretc|
|00001260| 68 65 64 20 76 61 6c 75 | 65 73 20 2a 2f 0d 0a 00 |hed valu|es */...|
|00001270| 00 00 00 09 70 74 72 61 | 6e 73 66 6f 72 6d 28 78 |....ptra|nsform(x|
|00001280| 2c 79 2c 64 78 2c 64 79 | 2c 74 72 61 6e 29 3b 0d |,y,dx,dy|,tran);.|
|00001290| 0a 00 00 00 00 7d 0d 0a | 00 00 00 00 0d 0a 00 00 |.....}..|........|
|000012a0| 00 00 2f 2a 20 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b |../* +++|++++++++|
|000012b0| 2b 2b 2b 2b 2b 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b |++++++++|++++++++|
|000012c0| 2b 2b 2b 2b 2b 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b |++++++++|++++++++|
|000012d0| 2b 2b 2b 2b 2b 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b |++++++++|++++++++|
|000012e0| 2b 2b 2b 2b 2b 20 2a 2f | 0d 0a 00 00 00 00 56 4f |+++++ */|......VO|
|000012f0| 49 44 20 70 72 69 6e 74 | 5f 68 69 73 74 6f 67 72 |ID print|_histogr|
|00001300| 61 6d 28 68 29 09 09 2f | 2a 20 50 72 69 6e 74 20 |am(h)../|* Print |
|00001310| 68 69 73 74 6f 67 72 61 | 6d 20 68 20 2a 2f 0d 0a |histogra|m h */..|
|00001320| 00 00 00 00 6c 6f 6e 67 | 20 2a 68 3b 0d 0a 00 00 |....long| *h;....|
|00001330| 00 00 7b 0d 0a 00 00 00 | 00 09 72 65 67 69 73 74 |..{.....|..regist|
|00001340| 65 72 20 69 6e 74 20 69 | 2c 63 6f 75 6e 74 3b 0d |er int i|,count;.|
|00001350| 0a 00 00 00 00 0d 0a 00 | 00 00 00 09 70 72 69 6e |........|....prin|
|00001360| 74 66 28 22 48 69 73 74 | 6f 67 72 61 6d 20 76 61 |tf("Hist|ogram va|
|00001370| 6c 75 65 73 3a 5c 6e 22 | 29 3b 0d 0a 00 00 00 00 |lues:\n"|);......|
|00001380| 09 63 6f 75 6e 74 20 3d | 20 30 3b 0d 0a 00 00 00 |.count =| 0;.....|
|00001390| 00 09 66 6f 72 20 28 69 | 20 3d 20 30 20 3b 20 69 |..for (i| = 0 ; i|
|000013a0| 20 3c 20 50 49 58 45 4c | 5f 53 49 5a 45 20 3b 20 | < PIXEL|_SIZE ; |
|000013b0| 69 2b 2b 29 20 7b 0d 0a | 00 00 00 00 09 09 70 72 |i++) {..|......pr|
|000013c0| 69 6e 74 66 28 22 25 36 | 6c 64 20 22 2c 2a 68 2b |intf("%6|ld ",*h+|
|000013d0| 2b 29 3b 0d 0a 00 00 00 | 00 09 09 69 66 20 28 21 |+);.....|...if (!|
|000013e0| 28 28 2b 2b 63 6f 75 6e | 74 29 25 38 29 29 20 70 |((++coun|t)%8)) p|
|000013f0| 72 69 6e 74 66 28 22 5c | 6e 22 29 3b 0d 0a 00 00 |rintf("\|n");....|
|00001400| 00 00 09 7d 0d 0a 00 00 | 00 00 7d 0d 0a 00 00 00 |...}....|..}.....|
|00001410| 00 0d 0a 00 00 00 00 0d | 0a 00 00 00 00 09 0d 0a |........|........|
|00001420| 00 00 00 00 2f 2a 20 3d | 3d 3d 3d 3d 3d 3d 3d 20 |..../* =|======= |
|00001430| 53 75 62 72 6f 75 74 69 | 6e 65 73 20 66 6f 72 20 |Subrouti|nes for |
|00001440| 73 65 74 74 69 6e 67 20 | 75 70 20 4c 55 54 53 2c |setting |up LUTS,|
|00001450| 20 69 66 20 79 6f 75 20 | 68 61 76 65 20 74 68 65 | if you |have the|
|00001460| 6d 20 3d 3d 3d 3d 3d 3d | 3d 3d 3d 20 2a 2f 0d 0a |m ======|=== */..|
|00001470| 00 00 00 00 0d 0a 00 00 | 00 00 56 4f 49 44 20 6c |........|..VOID l|
|00001480| 69 6e 5f 6c 75 74 73 28 | 29 09 09 2f 2a 20 53 65 |in_luts(|)../* Se|
|00001490| 74 20 4c 55 54 53 20 74 | 6f 20 6c 69 6e 65 61 72 |t LUTS t|o linear|
|000014a0| 20 28 67 72 65 79 2d 73 | 63 61 6c 65 29 20 2a 2f | (grey-s|cale) */|
|000014b0| 0d 0a 00 00 00 00 7b 0d | 0a 00 00 00 00 23 69 66 |......{.|.....#if|
|000014c0| 64 65 66 20 4c 55 54 53 | 0d 0a 00 00 00 00 09 72 |def LUTS|.......r|
|000014d0| 65 67 69 73 74 65 72 20 | 69 6e 74 20 69 3b 0d 0a |egister |int i;..|
|000014e0| 00 00 00 00 0d 0a 00 00 | 00 00 09 70 72 69 6e 74 |........|...print|
|000014f0| 66 28 22 2d 2d 20 4c 69 | 6e 69 61 72 69 7a 65 20 |f("-- Li|niarize |
|00001500| 6f 75 74 70 75 74 20 4c | 55 54 53 20 2d 2d 5c 6e |output L|UTS --\n|
|00001510| 22 29 3b 0d 0a 00 00 00 | 00 09 66 6f 72 20 28 69 |");.....|..for (i|
|00001520| 20 3d 20 30 20 3b 20 69 | 20 3c 20 50 49 58 45 4c | = 0 ; i| < PIXEL|
|00001530| 5f 53 49 5a 45 20 3b 20 | 69 2b 2b 29 20 7b 0d 0a |_SIZE ; |i++) {..|
|00001540| 00 00 00 00 09 09 77 72 | 69 74 65 5f 4c 55 54 28 |......wr|ite_LUT(|
|00001550| 52 45 44 2c 69 2c 28 50 | 49 58 45 4c 29 69 29 3b |RED,i,(P|IXEL)i);|
|00001560| 0d 0a 00 00 00 00 09 09 | 77 72 69 74 65 5f 4c 55 |........|write_LU|
|00001570| 54 28 47 52 45 45 4e 2c | 69 2c 28 50 49 58 45 4c |T(GREEN,|i,(PIXEL|
|00001580| 29 69 29 3b 0d 0a 00 00 | 00 00 09 09 77 72 69 74 |)i);....|....writ|
|00001590| 65 5f 4c 55 54 28 42 4c | 55 45 2c 69 2c 28 50 49 |e_LUT(BL|UE,i,(PI|
|000015a0| 58 45 4c 29 69 29 3b 0d | 0a 00 00 00 00 09 7d 0d |XEL)i);.|......}.|
|000015b0| 0a 00 00 00 00 23 65 6e | 64 69 66 0d 0a 00 00 00 |.....#en|dif.....|
|000015c0| 00 7d 0d 0a 00 00 00 00 | 0d 0a 00 00 00 00 2f 2a |.}......|....../*|
|000015d0| 20 2b 2b 2b 2b 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b | +++++++|++++++++|
|000015e0| 2b 2b 2b 2b 2b 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b |++++++++|++++++++|
|000015f0| 2b 2b 2b 2b 2b 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b |++++++++|++++++++|
|00001600| 2b 2b 2b 2b 2b 2b 2b 2b | 2b 2b 2b 2b 2b 2b 2b 2b |++++++++|++++++++|
|00001610| 2b 20 2a 2f 0d 0a 00 00 | 00 00 56 4f 49 44 20 70 |+ */....|..VOID p|
|00001620| 5f 73 70 65 63 74 72 75 | 6d 28 29 09 2f 2a 20 53 |_spectru|m()./* S|
|00001630| 65 74 20 4c 55 54 53 20 | 74 6f 20 61 20 73 70 65 |et LUTS |to a spe|
|00001640| 63 74 72 75 6d 20 28 70 | 73 65 75 64 6f 2d 63 6f |ctrum (p|seudo-co|
|00001650| 6c 6f 72 29 20 2a 2f 0d | 0a 00 00 00 00 7b 0d 0a |lor) */.|.....{..|
|00001660| 00 00 00 00 23 69 66 64 | 65 66 20 4c 55 54 53 0d |....#ifd|ef LUTS.|
|00001670| 0a 00 00 00 00 09 72 65 | 67 69 73 74 65 72 20 69 |......re|gister i|
|00001680| 6e 74 20 69 3b 0d 0a 00 | 00 00 00 09 50 49 58 45 |nt i;...|....PIXE|
|00001690| 4c 20 6a 2c 6b 3b 0d 0a | 00 00 00 00 09 69 6e 74 |L j,k;..|.....int|
|000016a0| 20 71 31 2c 71 32 2c 71 | 33 2c 71 34 3b 09 2f 2a | q1,q2,q|3,q4;./*|
|000016b0| 20 51 75 61 72 74 69 6c | 65 73 20 6f 66 20 4c 55 | Quartil|es of LU|
|000016c0| 54 20 72 61 6e 67 65 20 | 2a 2f 0d 0a 00 00 00 00 |T range |*/......|
|000016d0| 0d 0a 00 00 00 00 09 70 | 72 69 6e 74 66 28 22 2d |.......p|rintf("-|
|000016e0| 2d 20 53 65 74 20 6f 75 | 74 70 75 74 20 4c 55 54 |- Set ou|tput LUT|
|000016f0| 73 20 74 6f 20 61 20 73 | 70 65 63 74 72 75 6d 20 |s to a s|pectrum |
|00001700| 2d 2d 5c 6e 22 29 3b 0d | 0a 00 00 00 00 0d 0a 00 |--\n");.|........|
|00001710| 00 00 00 09 71 31 20 3d | 20 50 49 58 45 4c 5f 53 |....q1 =| PIXEL_S|
|00001720| 49 5a 45 3e 3e 32 3b 09 | 2f 2a 20 44 69 76 69 64 |IZE>>2;.|/* Divid|
|00001730| 65 20 4c 55 54 20 72 61 | 6e 67 65 20 62 79 20 34 |e LUT ra|nge by 4|
|00001740| 20 2a 2f 0d 0a 00 00 00 | 00 09 71 32 20 3d 20 71 | */.....|..q2 = q|
|00001750| 31 20 2b 20 71 31 3b 09 | 09 2f 2a 20 53 65 74 20 |1 + q1;.|./* Set |
|00001760| 75 70 20 71 75 61 72 74 | 69 6c 65 73 20 2a 2f 0d |up quart|iles */.|
|00001770| 0a 00 00 00 00 09 71 33 | 20 3d 20 71 32 20 2b 20 |......q3| = q2 + |
|00001780| 71 31 3b 0d 0a 00 00 00 | 00 09 71 34 20 3d 20 50 |q1;.....|..q4 = P|
|00001790| 49 58 45 4c 5f 53 49 5a | 45 3b 0d 0a 00 00 00 00 |IXEL_SIZ|E;......|
|000017a0| 0d 0a 00 00 00 00 2f 2a | 20 46 69 72 73 74 20 71 |....../*| First q|
|000017b0| 75 61 72 74 69 6c 65 20 | 6f 66 20 4c 55 54 20 72 |uartile |of LUT r|
|000017c0| 61 6e 67 65 2e 20 20 52 | 65 64 20 69 6e 63 72 65 |ange.  R|ed incre|
|000017d0| 61 73 65 73 20 6f 6e 6c | 79 20 2a 2f 0d 0a 00 00 |ases onl|y */....|
|000017e0| 00 00 09 66 6f 72 20 28 | 69 20 3d 20 30 20 3b 20 |...for (|i = 0 ; |
|000017f0| 69 20 3c 20 71 31 20 3b | 20 69 2b 2b 29 20 7b 0d |i < q1 ;| i++) {.|
|00001800| 0a 00 00 00 00 09 09 77 | 72 69 74 65 5f 4c 55 54 |.......w|rite_LUT|
|00001810| 28 52 45 44 2c 69 2c 28 | 50 49 58 45 4c 29 28 69 |(RED,i,(|PIXEL)(i|
|00001820| 3c 3c 32 29 29 3b 09 2f | 2a 20 52 65 64 20 72 61 |<<2));./|* Red ra|
|00001830| 6d 70 20 75 70 20 2a 2f | 0d 0a 00 00 00 00 09 09 |mp up */|........|
|00001840| 77 72 69 74 65 5f 4c 55 | 54 28 47 52 45 45 4e 2c |write_LU|T(GREEN,|
|00001850| 69 2c 4d 49 4e 50 49 58 | 29 3b 09 2f 2a 20 30 20 |i,MINPIX|);./* 0 |
|00001860| 67 72 65 65 6e 20 2a 2f | 0d 0a 00 00 00 00 09 09 |green */|........|
|00001870| 77 72 69 74 65 5f 4c 55 | 54 28 42 4c 55 45 2c 69 |write_LU|T(BLUE,i|
|00001880| 2c 4d 49 4e 50 49 58 29 | 3b 09 2f 2a 20 30 20 62 |,MINPIX)|;./* 0 b|
|00001890| 6c 75 65 20 2a 2f 0d 0a | 00 00 00 00 09 7d 0d 0a |lue */..|.....}..|
|000018a0| 00 00 00 00 2f 2a 20 53 | 65 63 6f 6e 64 20 71 75 |..../* S|econd qu|
|000018b0| 61 72 74 69 6c 65 20 6f | 66 20 4c 55 54 20 72 61 |artile o|f LUT ra|
|000018c0| 6e 67 65 2e 20 20 52 65 | 64 20 64 65 63 72 65 61 |nge.  Re|d decrea|
|000018d0| 73 65 73 2c 20 67 72 65 | 65 6e 20 69 6e 63 72 65 |ses, gre|en incre|
|000018e0| 61 73 65 73 20 2a 2f 0d | 0a 00 00 00 00 09 6a 20 |ases */.|......j |
|000018f0| 3d 20 4d 49 4e 50 49 58 | 3b 20 6b 20 3d 20 4d 41 |= MINPIX|; k = MA|
|00001900| 58 50 49 58 3b 0d 0a 00 | 00 00 00 09 66 6f 72 20 |XPIX;...|....for |
|00001910| 28 69 20 3d 20 71 31 20 | 3b 20 69 20 3c 20 71 32 |(i = q1 |; i < q2|
|00001920| 20 3b 20 69 2b 2b 29 20 | 7b 0d 0a 00 00 00 00 09 | ; i++) |{.......|
|00001930| 09 77 72 69 74 65 5f 4c | 55 54 28 52 45 44 2c 69 |.write_L|UT(RED,i|
|00001940| 2c 6b 29 3b 09 09 2f 2a | 20 52 65 64 20 64 65 63 |,k);../*| Red dec|
|00001950| 72 65 61 73 65 73 20 2a | 2f 0d 0a 00 00 00 00 09 |reases *|/.......|
|00001960| 09 6b 20 2d 3d 20 34 3b | 0d 0a 00 00 00 00 09 09 |.k -= 4;|........|
|00001970| 77 72 69 74 65 5f 4c 55 | 54 28 47 52 45 45 4e 2c |write_LU|T(GREEN,|
|00001980| 69 2c 6a 29 3b 09 09 2f | 2a 20 47 72 65 65 6e 20 |i,j);../|* Green |
|00001990| 69 6e 63 72 61 73 65 73 | 20 2a 2f 0d 0a 00 00 00 |incrases| */.....|
|000019a0| 00 09 09 6a 20 2b 3d 20 | 34 3b 0d 0a 00 00 00 00 |...j += |4;......|
|000019b0| 09 09 77 72 69 74 65 5f | 4c 55 54 28 42 4c 55 45 |..write_|LUT(BLUE|
|000019c0| 2c 69 2c 4d 49 4e 50 49 | 58 29 3b 09 2f 2a 20 30 |,i,MINPI|X);./* 0|
|000019d0| 20 62 6c 75 65 20 2a 2f | 0d 0a 00 00 00 00 09 7d | blue */|.......}|
|000019e0| 0d 0a 00 00 00 00 2f 2a | 20 54 68 69 72 64 20 71 |....../*| Third q|
|000019f0| 75 61 72 74 69 6c 65 20 | 6f 66 20 4c 55 54 20 72 |uartile |of LUT r|
|00001a00| 61 6e 67 65 2e 20 20 47 | 72 65 65 6e 20 64 65 63 |ange.  G|reen dec|
|00001a10| 72 65 61 73 65 73 2c 20 | 42 6c 75 65 20 69 6e 63 |reases, |Blue inc|
|00001a20| 72 65 61 73 65 73 20 2a | 2f 0d 0a 00 00 00 00 09 |reases *|/.......|
|00001a30| 6a 20 3d 20 4d 49 4e 50 | 49 58 3b 20 6b 20 3d 20 |j = MINP|IX; k = |
|00001a40| 4d 41 58 50 49 58 3b 0d | 0a 00 00 00 00 09 66 6f |MAXPIX;.|......fo|
|00001a50| 72 20 28 69 20 3d 20 71 | 32 20 3b 20 69 20 3c 20 |r (i = q|2 ; i < |
|00001a60| 71 33 20 3b 20 69 2b 2b | 29 20 7b 0d 0a 00 00 00 |q3 ; i++|) {.....|
|00001a70| 00 09 09 77 72 69 74 65 | 5f 4c 55 54 28 52 45 44 |...write|_LUT(RED|
|00001a80| 2c 69 2c 4d 49 4e 50 49 | 58 29 3b 09 2f 2a 20 30 |,i,MINPI|X);./* 0|
|00001a90| 20 72 65 64 20 2a 2f 0d | 0a 00 00 00 00 09 09 77 | red */.|.......w|
|00001aa0| 72 69 74 65 5f 4c 55 54 | 28 47 52 45 45 4e 2c 69 |rite_LUT|(GREEN,i|
|00001ab0| 2c 6b 29 3b 09 09 2f 2a | 20 47 72 65 65 6e 20 64 |,k);../*| Green d|
|00001ac0| 65 63 72 65 61 73 65 73 | 20 2a 2f 0d 0a 00 00 00 |ecreases| */.....|
|00001ad0| 00 09 09 6b 20 2d 3d 20 | 34 3b 0d 0a 00 00 00 00 |...k -= |4;......|
|00001ae0| 09 09 77 72 69 74 65 5f | 4c 55 54 28 42 4c 55 45 |..write_|LUT(BLUE|
|00001af0| 2c 69 2c 6a 29 3b 09 09 | 2f 2a 20 42 6c 75 65 20 |,i,j);..|/* Blue |
|00001b00| 69 6e 63 72 65 61 73 65 | 73 20 2a 2f 0d 0a 00 00 |increase|s */....|
|00001b10| 00 00 09 09 6a 20 2b 3d | 20 34 3b 0d 0a 00 00 00 |....j +=| 4;.....|
|00001b20| 00 09 7d 0d 0a 00 00 00 | 00 2f 2a 20 46 6f 72 74 |..}.....|./* Fort|
|00001b30| 68 20 71 75 61 72 74 69 | 6c 65 20 6f 66 20 4c 55 |h quarti|le of LU|
|00001b40| 54 20 72 61 6e 67 65 2e | 20 20 42 6c 75 65 20 64 |T range.|  Blue d|
|00001b50| 65 63 72 65 61 73 65 73 | 20 2a 2f 0d 0a 00 00 00 |ecreases| */.....|
|00001b60| 00 09 6b 20 3d 20 4d 41 | 58 50 49 58 3b 0d 0a 00 |..k = MA|XPIX;...|
|00001b70| 00 00 00 09 66 6f 72 20 | 28 69 20 3d 20 71 33 20 |....for |(i = q3 |
|00001b80| 3b 20 69 20 3c 20 71 34 | 20 3b 20 69 2b 2b 29 20 |; i < q4| ; i++) |
|00001b90| 7b 0d 0a 00 00 00 00 09 | 09 77 72 69 74 65 5f 4c |{.......|.write_L|
|00001ba0| 55 54 28 52 45 44 2c 69 | 2c 4d 49 4e 50 49 58 29 |UT(RED,i|,MINPIX)|
|00001bb0| 3b 09 2f 2a 20 30 20 72 | 65 64 20 2a 2f 0d 0a 00 |;./* 0 r|ed */...|
|00001bc0| 00 00 00 09 09 77 72 69 | 74 65 5f 4c 55 54 28 47 |.....wri|te_LUT(G|
|00001bd0| 52 45 45 4e 2c 69 2c 4d | 49 4e 50 49 58 29 3b 09 |REEN,i,M|INPIX);.|
|00001be0| 2f 2a 20 30 20 67 72 65 | 65 6e 20 2a 2f 0d 0a 00 |/* 0 gre|en */...|
|00001bf0| 00 00 00 09 09 77 72 69 | 74 65 5f 4c 55 54 28 42 |.....wri|te_LUT(B|
|00001c00| 4c 55 45 2c 69 2c 6b 29 | 3b 0d 0a 00 00 00 00 09 |LUE,i,k)|;.......|
|00001c10| 09 6b 20 2d 3d 20 34 3b | 0d 0a 00 00 00 00 09 7d |.k -= 4;|.......}|
|00001c20| 0d 0a 00 00 00 00 23 65 | 6e 64 69 66 0d 0a 00 00 |......#e|ndif....|
|00001c30| 00 00 7d 0d 0a 00 00 00 | 00 0d 0a 00 00 00 00 2f |..}.....|......./|
|00001c40| 2a 20 3d 3d 3d 3d 3d 3d | 3d 3d 3d 3d 3d 3d 3d 3d |* ======|========|
|00001c50| 3d 3d 20 45 6e 64 20 6f | 66 20 73 69 6d 73 75 62 |== End o|f simsub|
|00001c60| 73 2e 63 20 3d 3d 3d 3d | 3d 3d 3d 3d 3d 3d 3d 3d |s.c ====|========|
|00001c70| 3d 3d 3d 3d 20 2a 2f 0d | 0a 00 00 00 00 0d 0a 00 |==== */.|........|
|00001c80| 00 00 00 2f 2a 20 3c 2d | 2d 20 46 49 4c 45 20 42 |.../* <-|- FILE B|
|00001c90| 52 45 41 4b 20 2d 2d 3e | 20 2a 2f 0d 0a 00 00 00 |REAK -->| */.....|
|00001ca0| 00 25 20 1a 00 00 23 65 | 6e 64 69 66 0d 0a 00 00 |.% ...#e|ndif....|
|00001cb0| 00 00 7d 0d 0a 00 00 00 | 00 0d 0a 00 00 00 00 2f |..}.....|......./|
|00001cc0| 2a 20 3d 3d 3d 3d 3d 3d | 3d 3d 3d 3d 3d 3d 3d 3d |* ======|========|
|00001cd0| 3d 3d 20 45 6e 64 20 6f | 66 20 73 69 6d 73 75 62 |== End o|f simsub|
|00001ce0| 73 2e 63 20 3d 3d 3d 3d | 3d 3d 3d 3d 3d 3d 3d 3d |s.c ====|========|
|00001cf0| 3d 3d 3d 3d 20 2a 2f 0d | 0a 00 00 00 00 0d 0a 00 |==== */.|........|
+--------+-------------------------+-------------------------+--------+--------+